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We have performed a simulation study of Langmuir turbulence in the Earth's ionosphere by means 
of a Zakharov model with parameters relevant for the F layer. The model includes dissipative terms 
to model collisions and Landau damping of the electrons and ions, and a linear density profile, which 
models the ionospheric plasma inhomogeneity whose length scale is of the order 10-100 km. The 
injection of energy into the system is modeled by a constant source term in the Zakharov equation. 
Langmuir turbulence is excited "on-demand" in controlled ionospheric modification experiments 
where the energy is provided by an HF radio beam injected into the overhead ionospheric plasma. 
The ensuing turbulence can be studied with radars and in the form of secondary radiation recorded 
by ground-based receivers. We have analyzed spectral signatures of the turbulence for different 
sets of parameters and different altitudes relative to the turning point of the linear Langmuir mode 
where the Langmuir frequency equals the local plasma frequency. By a parametric analysis, we have 
derived a simple scaling law, which links the spectral width of the turbulent frequency spectrum to 
the physical parameters in the ionosphere. The scaling law provides a quantitative relation between 
the physical parameters (temperatures, electron number density, ionospheric length scale, etc.) and 
the observed frequency spectrum. This law may be useful for interpreting experimental results. 



I. INTRODUCTION 

The interaction between an electromagnetic wave and 
an inhomogeneous plasma layer leads to phenomena on 
differe nt temporal and sp atial scales. Laboratory experi- 
ments [Kim et al\ , Il974| have demonstrated cavity for- 
mation and trapping of radio frequency (RF) electro- 
static fields due to density modification of the plasma 
near the critical density where the plasma frequency 
equals that of the electrostatic RF pump field. Iono- 
spheric cavitons were observed near the critical layer of 
a high- power radio wave in jected into the overhead iono- 
sphere [Wong et all Il987j . In these experiments, both 
decimeter-scale cavitons due to the ponderomotive force 
by mode-converted electrostatic waves (electrostatic cavi- 
tons) and kilometer-sized cavitons due to the thermal 
pressure of the heated electrons (thermal cavitons) were 
observed. 

By analyzing the frequency spectrum of the electro- 
magnetic waves recorded by ground-based receivers, it 
was found experimentally in Troms0 that strong, system- 
atic, structured, wide-band secondary HF radiation es- 
capes from the turbulent interaction region Thide et~al\ , 
Il982| |. This and other observations demonstrated that 
complex interactions, incl uding weak and str o ng EM 
turbu lence, are involved [Carozzi et al. |, 120021: \Levsel 
l200l| . Naturally enhanced ion acou stic lines (NEIALs 
observed by several radar facilities Foster et all Il98 
\Sedaemore-Schulthess and St. Maurice] , 200l| have in 
some cases been interp r eted as the result of Lang- 
muir turbulence \Forme\, 1993t Guio and Formi . [2006; 



simulation studies in homogeneous plasmas have demon- 
strated the ba sic processes of Langmuir t urbulence and 
wave collapse WuBois et all 1991 . 1993 : Hansen et 



1992t \Mi0lhulWoO& \Nicholson et all 11984 \Robins~ . 

Simulation studies with an inhomogeneous Za- 
kharov equation, where the electromagnetic field was 
modeled by an external dipolar field, demonstrated den- 
sity modification and the excitation of localized clcc- 
tric fields and ion-acou s tic os cillations in the plasma 
[Morales and Lei . 11974 Il977 |. Resonant absorption 



I Sedaemore-Schulthess and St. Ma.unrA b r^S^al 



of electromagnetic waves and electrostatic turbulence 
was s tudied experimentally and numerically [^2S^eFjzI|, 
Il991 | with a modified Zakharov model and demonstrated 
the density modification and transition to turbulence for 
different amplitudes of the electromagnetic RF pump 



In the present article, we present results from long-time 
series simulations of the driven Zakharov system of equa- 
tions in order to investigate the statistical properties of 
Langmuir turbulence in an inhomogeneous plasma, and 
to obtain the dependence of the spectrum on physical 
parameters. The manuscript is organized in the follow- 
ing fashion. In Section 2, we present the mathematical 
model in the form of the driven Zakharov equations for 
an inhomogeneous plasma background, including Landau 
damping and collision terms. The numerical setup and 
results are presented in Section 3, where frequency spec- 
tra and their dependence on physical variables are dis- 
cussed, and an approximate scaling law for the spectral 
width is derived. The results obtained are summarized 
in Section 4. 
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II. THE INHOMOGENEOUS ZAKHAROV 
EQUATION 

In the Zakharov model, the Langmuir wave field in the 
presence of a linear density gradient, a constant dipolc 
pump fiel d and slowly varying density fl uctuations is gov- 
erned by [Morales and -Led . 11974 119771 ] 
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where i is the imaginary unit, io pe — (n^e 2 / e^mf,) 1 ' 2 
is the electron plasma frequency, x is the vertical spa- 
tial coordinate, L is the local length scale of the iono- 
spheric profile, n s is the slowly varying electron density 
perturbation due to ion-acoustic fluctuations, no is the 
equilibrium electron number density, 7 l is the collision 
frequency, \oe — VTe/^pe is the electron Debye length, 
VTe = (kBTe/rrie) 1 / 2 is the electron thermal speed, e is 
the magnitude of the electron charge, m e is the electron 
mass, Eg is the vacuum permittivity, kg is Boltzmann's 
constant, and T e is the electron temperature. The source 
term -Epump repr esents the energy provided by electro- 
magnetic waves [Thide et all Il982j . If the radio wave 
has an oblique incidence to the plasma gradient, then it 
will be reflected at an altitude slightly lower than the crit- 
ical height where the pump frequency equals the plasma 
frequency. A small fraction of the electromagnetic wave 
will tunnel up to the critical height where it excites 
large-amplitude electrostatic waves. In our model, the 
dipole pump represents the evanescent electromagnetic 
wave that has tunneled up to the critical altitude. Other 
mode conversion processes may b e important du e to the 
presence of the geomagnetic field Mj0lhui , Il99fjj | , which 
we have neglected here. The ion acoustic fluctuations in 
the presence of the ponderomotive force are governed by 
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FIG. 1: a) Profiles of the amplitude of Langmuir electric field, 
\E\ and b) the electron number density perturbation n s /no 



for E pun 



1.0 x 10 V/m. c) and d): The same as a) and 



b) but for -Epump = 1.0 x 10 3 V/m. The other parameters 
are L = 50 km, -y L = 1.0 x 10 3 s~\ and 7s = 1.0 x 10 3 s _1 . 
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FIG. 2: The time-dependent amplitude at z — of the Lang- 
muir electric field \E\ (upper panel) and of the electron num- 
ber density perturbation n s /no (lower panel) for E pn mp = 
1.0 x 10 _4 V/m. The other parameters are L = 50 km, 
7l = 1.0 x 10 3 s"\ and 7„ = 1.0 x 10 3 s" 1 . 



where C s — [fcs(T e + ST^/nii] 1 / 2 is the ion acoustic 
speed, Tj is the ion temperature, m, is the ion mass, and 
7 S is the ion collision frequency Nicholson et all 11984 
due to collisions and/or Landau damping, which, for sim- 
plicity, we have taken to be constant. 



III. NUMERICAL RESULTS 

In order to investigate the nonlinear dynamics and 
spectral features of the Zakharov system, we have solved 
Eqs. |U) and © for different sets of parameters. We 
have used n = 5 x KFnT 3 , T e = T, = 2000 K, and 
we have used n%i — 26.8 x 10~ 27 kg (oxygen ions) in 

VTe = 



our simulations, so that uj pe = 4.0 x 10 s 
1.7 x 10 5 m/s and C s = 2.0 x 10 3 m/s. The pump 
frequency, which is set equal to the local plasma fre- 
quency, is fo = (o pe /2iT = 6.34 MHz. In the numerical 



procedure, the Zakharov system was advanced in time 
with a standard fourth-order Runge-Kutta scheme with a 
timestep of 10 -6 s. The spatial derivatives were approxi- 
mated with centered, second-order difference approxima- 
tions with a spatial grid spacing of 5 cm. For weaker 
pump fields Ep Ump of order 10~ 4 V/m, we used the spa- 
tial domain —100 m < z < 100 m, and for the stronger 
pump fields of order 10~ 3 V/m, we used the spatial do- 
main -200 m < z < 200 m. 

In Fig. [TJ we show the profiles of the electric field and 
the density fluctuations at the end of the simulation, for 
two different values of the pump field. For a pump field of 
-Epump = O.lmV/s [panels a) and b)], we see an excited 
electrostatic field of the order \E\ ~ 3 V/m correlated 
with a density depletion of n s /no ~ —4 x 10~ 4 . The 
turbulent field is localized in a relatively small interval 
of z between —10 m and 10 m, correlated with electric 
field maxima are density minima. For the larger pump 
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FIG. 3: Frequency spectra (in dB) of E at different altitudes, 
a) z — 10 m, z = 5 m, z — m, z = —5 m, and z = —10 m, 
(top to bottom panel), for E pulnp — 10 _4 V/m and b) z = 
30m, 2 = 15m, z = 0m, z = —25m, and z = —50m, (top to 
bottom panel), for -E pump = 10~ 3 V/m. The other parameters 
are L = 50 km, j L = 1.0 x 10 3 s~\ and j s = 1.0 x 10 3 s _1 . 



field -Epump = l.OV/m, displayed in panels c) and d), 
the turbulence takes place in a larger altitude interval, 
between z ~ —60 m and z ~ 30 m. Here, the Langmuir 
field is larger, of the order 10-20 V/m, and the maximum 
density depletion, correlated with electric field maxima, 
is of the order n s /n ~ — 10~ 2 . 

The time dependent amplitude of the Langmuir field 
and the density fluctuations at z = are shown in Fig. [21 
for the first 50 ms (physical time) . After an initial growth 
phase of the fields during the first 3 ms, the fluctuations 
enter a steady-state turbulent phase during the rest of 
the time. In Fig. [31 we are analyzing the frequency spec- 
trum of the electric field at different altitudes. In order to 
obtain the frequency spectrum at each altitude, the sim- 
ulation was run for 500 ms. The time series of the elec- 
tric field was subdivided into five time-slices of 100 ms, 
and each time slice is Fourier transformed in time with 
a Hamming window. Finally, the average is taken of the 
five power spectra. The averaged power spectrum (in 
dB) is visualized in Fig. [3] for E pU mp = O.lmV/m and 
Epump = 1 mV/m, at different altitudes. We see that es- 
pecially the the down-shifted part of the power spectrum 
differs significantly between different altitudes, while the 
upshifted part of the spectra vary in their amplitudes but 
their logarithmic slopes are almost the same at different 
altitudes. In column a) of Fig. 3, we see that the spec- 
trum is much wider at z — 5 m than at z = —5 m. Hence, 
in situ measurements by satellites or rockets would be 
very sensitive to the exact location of the measuring 
probe. In a thought ground-based experiment, the elec- 
trostatic turbulence converts to electromagnetic radia- 



tion by some process [Stubbe et all 1 1984 1 , and the escap- 
ing radiation is monitored by receiving equipment of the 
ground. The conversion process may involve the mode- 
conversion of the electrostatic waves against density gra- 
dients, resulting in a mixture of frequency components 
from different altitudes in the resultant electromagnetic 
radiation. The radiation field from an extended body can 
be expressed in terms of an integral over the retarded 
current density. For electrostatic turbulence, the current 
density and electrostatic field are linked by the Maxwell 
equation d~E/dt = — j/eo, an d hence the radiation field 
is proportional to the spatial integral of the electrostatic 
field over the turbulent region. We here do not model in 
detail the conversion process but will assume that the es- 
caping radiation contains frequency components from the 
spatial average of the electrostatic field (possibly modi- 
fied by the propagation characteristics of the ionospheric 
plasma). This model may be valid as long as the wave- 
lengths of the electrostatic waves are much smaller than 
that of the vacuum electromagnetic waves. In our case 
the electrostatic waves are on decimeter or meter scale, 
while the electromagnetic waves are ~ 60 m. Hence, we 
have taken the spatial average of the complex electric 
field E over the simulation box. The spatially averaged 
electric field is Fourier analyzed in a similar manner as 
in Fig. [31 and the resulting frequency spectra are plotted 
in Figs. 0H7] for different sets of parameters. 

In principle, the eight physical parameters that can be 
varied are Epump, E, 7i, Is-, T e , T u n , and mj. However, 
the number of parameters can b e reduced by normalizin g 
the system, in the spirit of, e.g.. \Morales and Led |l977 |. 
Introducing dimensionless, primed variables according to 
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ls = (C 2 s co pe /3v 2 Te H, (10) 
the Zakharov system ((T])-(G3) is transformed into 



<9E' 



rP- F" 

pump V / 



4 



d 2 n' s ,dn' s 8 2 n' s d 2 \E'\ 



'12 



at' 2 



df 



dx' 2 



dx' 2 



(12) 



and we see that the scaled system have only four free 
parameters L', j' L , E pump , and j' s . We can state the 
similarity principle that for any change of parameters 
(temperatures, densities, etc.) in the dimensional system 
such that the dimensionless parameters V , j' L , E' pump , 
and 7g are unchanged, the system is self-similar, i.e., it 
will have the same solutions as the original system up 
to linear scalings. We have chosen to vary only the four 
parameters E pump , L, 7^ and 7 S in our numerical simu- 
lations, while keeping the other parameters constant. 
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FIG. 4: a) Frequency spectra (in dB) of E (spatially aver- 
aged) for different pump amplitudes -Ep Um p = 0.5 x 1CP 4 V/m, 
£ pump = 0.8 x 10" 4 V/m, S pump = 1.0 x 10~ 4 V/m, E pump = 
1.2 x 10~ 4 V/m, and E pump = 1.5 x 10~ 4 V/m (top to bot- 
tom panel), for L = 50km. The spectral width is denned 
as the frequency difference between the high-frequency and 
low-frequency points where the power has dropped 30 dB 
compared to the downshifted maximum, as indicated with 
dashed lines for the case E puulp — 1.0 x 10 -4 V/m; here 
the spectral width is approximately 5 kHz. b) Frequency 
spectra (in dB) of E (spatially averaged) for different iono- 
spheric length scales L = 15 km, L = 25 km, L = 50 km, 
L = 100km, and L = 150km (top to bottom panel), for 
-Epump = 1-0 x 10 _4 V/m. The other parameters are 71, = 
1.0 x 10 3 s" 1 and 7 S = 1.0 x 10 3 s" 1 

All the spectra in Figs. 0H7J have some common fea- 
tures. The spectra are asymmetric, where the down- 
shifted part of the frequency spectrum is about three 
times wider than the upshifted one, although the rela- 
tive width seems to increase for the downshifted spectrum 
for larger pump amplitude; see Figs. [S] and [7J We also 
see a downshifted maximum in the spectra, whose down- 
shift varies from a few kHz up to ~ 5 kHz, depending 
on the parameters. These general s pectral featu r es are 
frequently observed in experiments [Frolov et all l2004t 
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FIG. 5: a) Frequency spectra (in dB) of E (spatially averaged) 
for electron collision frequencies jl = 2.0 x 10 3 s _1 , jl = 
1.5 x 10 3 s"\ 7 L = 1.0 x 10 3 s" 1 , 7 L = 0.5 x 10 3 s -1 , and -y L = 
0.2 x 10 3 s~ 1 (top to bottom panel), for 7 S = 1.0 x 10 3 s _1 . 
b) Frequency spectra (in dB) of E (spatially averaged) for 
different ion collision frequencies 7l = 2.0 x 10 3 s _1 , 7s = 
1.5 x 10 3 s"\ 7 S = 1.0 x 10 3 s~\ 7 S = 0.5 x 10 3 s~\ and 7s = 
0.2 x 10 3 s _1 (top to bottom panel), for -jl 
The other parameters are L 
10" 4 V/m. 



50 km and E p 



1.0 x 



\Thide et all Il982l . Il989t I Thidi Il990j . A few examples of 
measured SEE spectra are shown in Fig. [8] below, where 
we see asymmetric spectra with downshifted components. 
We also observe a slow temporal change of the spectrum, 
which indicates that the plasma parameters in the iono- 
sphere are slowly modified by the large-amplitude elec- 
tromagnetic wave. 

In Fig.[U we have varied the pump field -E pU mp and the 
ionospheric length scale L. We define the spectral width 
as the frequency difference between the high-frequency 
and low-frequency points where the power has dropped 
30 dB compared to the downshifted maximum, as in- 
dicated with dashed lines in column a) of Fig. 2] for 
-Epump = 1-0 x 1CP 4 V/m; here the spectral width is 
approximately 5 kHz. We see in column a) of Fig.[2]that 
the spectral width is almost linearly proportional to the 
amplitude of E pump . The spectra also widens for larger 
ionospheric length scales, seen in column b) of Fig.[4j and 
the spectral width is approximately proportional to the 
square root of the length scale L. In Fig. [5J we have var- 
ied the collision frequencies 7^ and 7s- We here see that 
the width of the spectrum increases for smaller values of 
7l and 7 S , and is approximately inversely proportional 
to the square root of jl and 7 S . 

In Figs. [S] and [71 we have repeated the simulations for a 
one order of magnitude larger pump field than in Figs. 0] 
and0 The qualitative results are the same in Figs.[5]and 
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FIG. 6: a) Frequency spectra (in dB) of E (spatially aver- 
aged) for different pump amplitudes -Ep Ump = 0.5 x 10~ 3 V/m, 
Spump = 1.0 x 10~ 3 V/m, and E pump = 2 x 10~ 3 V/m (top 
to bottom panel), for L — 50km. We have indicated the 
spectral width with dashed lines, in the same manner as in 
Fig- IU for the case E pump = 1.0 x 10~ 3 V/m; here the spec- 
tral width is approximately 50 kHz. b) Frequency spectra (in 
dB) of E (spatially averaged) for different ionospheric length 
scales L = 25 km, L = 50 km, and L — 100 km (top to bottom 
panel), for E pump = 1.0 x 10~ 3 V/m. The other parameters 
are jt = 1.0 x 10 3 s~ 1 and 7a = 1.0 x 10 3 s~\ 
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FIG. 7: a) Frequency spectra (in dB) of E (spatially averaged) 
for electron collision frequencies jl = 2.0 x 10 3 s _1 , 7l = 



1.0 x 10 3 s _1 , and = 0.5 x 10 3 s _1 , (top to bottom panel), 
for 7a = 1.0 x 10 3 s _1 . b) Frequency spectra (in dB) of E 
(spatially averaged) for different ion collision frequencies 7 S = 
2.0xl0 3 s~\ 7s = l.OxlO 3 s _1 , and 7s = 0.5xl0 3 s _1 , (top to 
bottom panel), for 7l = 1.0 x 10 3 s _1 . The other parameters 
are L = 50 km and -E pump = 1.0 x 10~ 3 V/m. 



[Jj In Fig. [6] we see that the with of the spectrum is ap- 
proximately linearly proportional to E pump . Comparing 
the case E pump = 1CP 3 V/m in column a) of Fig.[6l where 
the spectral with is 50 kHz with the case -E pU mp = 10~ 4 
V/m in column a) of Fig. [4] where the spectral width is 
5 kHz, we see that the linear dependence of the spectral 
width on the pump amplitude holds over one decade of 
pump amplitudes. Similar to the results in Figs. 0] and 
02 the spectra in Figs. [6] and [7J widen for increasing L 
and decreasing 7^ and 7s , though the downshifted parts 
of the spectra are slightly less sensitive to variations in 
these parameters. Hence, the spectral width obeys the 
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FIG. 8: a) SEE spectra recorded at the Sura heating fa- 
cility in Nizhny Novgorod, Russia, on 25 September 1991 at 
13:30-13:50 LT. The spectra are measured at different times 
(indicated in the boxes) after the switch-on of the pump wave, 
in a 10 s on, 50 s off mode with pump frequency 5828 MHz 
and ERP power 30 MW. The spectra are primarily down- 
shifted, and a slow modification of the spectra can be seen, 
with wider spectra at later times. The narrow spikes a re due 
to interfering radio stations. (After {prolov et ~al\ [2004 ].) 



approximate scaling law 
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where C is a normalization constant that depends on 
temperatures, the density, etc., but not on -Epump, L, 7^ 
or 7 S . The scaling law can also be expressed in terms of 
the normalized, primed variables in Eq. (|12p as 



[A/'] 



E' L' 1 ' 2 

^,1 "pump-" 

11/2 11/2 - 
1l Is 



(14) 



where C is a numerical scaling factor that does not de- 
pend on any physical parameters in the Zakharov system. 

By the scaling of the time ([3|) in the normalized system 
(fT2)) . we have the the relation 



dv Te 
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between the dimensional and dimensionless spectral 
widths. Inserting (fT3")) and (jT5J) into (fT5")) , and elim- 
inating E ] 



L, j L and 7 S with the help of 0- 
(IT0|) . we find, after some algebra, the relation C = 



pump ? 
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[(\/6/18)^ 2 C s 5/2 e/4 e (™ e ^) 1/2 ]C" between the di- 
mensional and dimensionless scaling factors. Hence, the 
scaling law (flU|) can be expressed as 
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where A' = (%/6/18)C" is a numerical, dimensionless fac- 
tor, which is determined by the definition of the spectral 
width. As an example, we take the third row from top 
of panels in Fig. 2J where we see that the spectral width 
is [A/] w ~ 5 x 10 3 Hz, and where the parameters are 



10" 4 V/m, L = 50 x 10 3 m, 



1l = Is 



10"V 



etc. Inserting the spectral width and the parameters into 
Eq. (|16jl. we determine the scaling factor to A! = 4. We 
have thus in Eq. (| 16[) a quantitative measure of the de- 
pendence of the frequency spectrum on the physical pa- 
rameters in the ionosphere. 

Returning to the experimental results in Fig. [51 we see 
that the spectral width of the SEE spectrum increases 
slowly with time, from approximately 40 kHz in the time 
intervals during the first second, to approximately 80-90 
kHz after 2-4 seconds. This could possibly be due to 
an increase of the electron temperature, so that T e 3> 
Tj, and a resulting decrease of the ion Landau damping 
7s 



Krall and Trivelviec4 . Il973 | with time, which in turn 



would lead to an increase of the spectral width according 
to our scaling law (fTBjl . 



IV. SUMMARY 

In conclusion, we have performed a simulation study of 
Langmuir turbulence with parameters relevant for the F 
layer in the Earth's ionosphere. We have used a Zakharov 
model, including collisions and an ionospheric plasma ho- 
mogeneity, whose length scale is of the order 10-100 km. 
The energy sourc e is assumed to be provided by electro- 
magnetic waves [Thide et all Il982j in controlled iono- 
spheric wave injection experiments, and is modeled by a 
constant source term in the Zakharov equation. We have 
analyzed the frequency spectra for different sets of pa- 
rameters and different altitudes relative to the classical 
turning point of the Langmuir waves. By a paramet- 
ric study, we have derived a simple scaling law, which 
links the spectral width of the turbulent frequency spec- 
trum to the physical parameters in the ionosphere. The 
scaling law provides a quantitative relation between the 
physical parameters (temperatures, electron number den- 
sity, ionospheric length scale, etc.) and the observed fre- 
quency spectrum. An application of our results is to 
compare our scaling law with SEE spectra, with the as- 
sumption that the processes creating the SEE reflect the 
spectral structure of the turbulent E field via mode con- 
version of electrostatic waves to electromagnetic radia- 
tion. 
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